# read in data
d1 <- read.csv("multiTimeline.csv",skip = 2)
names(d1) <- c("week","ej")
# generate linear model
n = 1:length(d1$ej)
m1 = lm(d1$ej~n)
# plot data
p1 <- plot_ly(data = d1, x =~week, y =~ej, type = 'scatter', mode = 'lines+markers', name= "ej")%>%
add_lines(data = d1, x = ~week, y =predict(m1), type = 'line', name="Linear Model")
p1